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ЛЕТАТЕЛЬНОГО АППАРАТА 


(57) Реферат: 

Изобретение относится к авиационной технике. 
Способ заключается в том, что измеряют 
вертикальное ускорение и с помощью датчика 
статического давления приборную 
барометрическую высоту, вертикальную скорость 
получают интегрированием разности вертикального 
ускорения и поправки ускорения, барометрическую 
высоту - интегрированием разности вертикальной 
скорости и поправки скорости, разность между 
вычисленной барометрической ВЫСОТОЙ и 


приборной барометрической высотой используют 
для вычисления упомянутых поправок, при этом 
вертикальную скорость умножают на безразмерный 
множитель, равный отношению температуры 
воздуха по стандартной атмосфере, вычисленной 
по значению  барометрической высоты, к 
измеренной температуре воздуха. Изобретение 
позволяет повысить точность определения 
барометрической высоты и вертикальной скорости 
в условиях летательного аппарата в нестандартной 
атмосфере. 4 ил. 
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(54) METHOD OF DETERMINATION OF BAROMETRIC ALTITUDE AND VERTICAL VELOCITY OF 


FLYING VEHICLE 


(57) Abstract: 

FIELD: aeronautical engineering. 

SUBSTANCE: proposed method consists in 
measuring the vertical acceleration and indicated 
barometric altitude by means of static pressure 


sensor; vertical velocity is X obtained through 
integration of difference of vertical 
acceleration and acceleration correction; 
barometric altitude is obtained through 
integration of difference of vertical velocity 
and velocity correction; difference between 
calculated barometric altitude and indicated 


barometric altitude is used for calculation of 
said corrections; vertical velocity is multiplied 
by unitless multiplier equal to ratio of air 
temperature by standard atmosphere calculated by 
magnitude of barometric altitude to measured 


temperature of air. 

EFFECT: enhanced accuracy of determination of 
barometric altitude and vertical velocity under 
off-standard conditions of flight. 

4 dwg 


Фиг. 1 (прототип) 
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Изобретение относится к авиационной технике и может быть использовано для 
повышения точности определения на борту летательного аппарата вертикальной скорости 
и барометрической высоты при полете в условиях нестандартной атмосферы. 

Используемые сокращения и термины: 

ЛА - летательный аппарат, 

СВС - система воздушных сигналов, 

ИНС - инерциальная система. 

Прототипом изобретения является способ определения барометрической высоты и 
вертикальной скорости, на котором основано функционирование устройства, описанного в 
патенте США №4882697 (кл. 364/454) в качестве уровня техники. Он использует 
комплексирование сигнала вертикального ускорения от ИНС и сигнала приборной 
барометрической высоты, вычисляемой по сигналу датчика статического давления. Для 
обработки упомянутых сигналов применяют алгоритм «баро-инерциальный фильтр», 
структурная схема которого приведена на фигуре 1. На схеме приняты следующие 
обозначения: 

ось Z - ось, направленная вертикально вверх, 

AZ, BZ - входные сигналы вертикального ускорения и приборной барометрической 
высоты, 

VZ, PZ - выходные сигналы оценок вертикальной скорости и барометрической высоты 
фильтра, 

1 - интеграторы, ® - сумматоры, 
g 


G1, G2 n G3 - коэффициенты усиления. 

C помощью ИНС измеряют вертикальное ускорение, а с помощью датчика статического 
давления - барометрическую высоту. Выходной сигнал вертикальной скорости VZ 
формируют интегрированием вертикального ускорения. При этом вводят коррекцию и в 
качестве интегрируемой величины используют разность между сигналом вертикального 
ускорения AZ и сигналом коррекции ускорения $7. Выходной сигнал высоты PZ формируют 
интегрированием разности сигнала вертикальной скорости \ и сигнала коррекции 
вертикальной скорости S5. Сигналы коррекции ускорения S7 и вертикальной скорости S5 
формируют по сигналу рассогласования Е, представляющего собой разность между 
сигналами S6 (он же PZ) и BZ. Составляющую $3 сигнала коррекции ускорения 57 и сигнал 
коррекции вертикальной скорости $5 формируют умножением сигнала рассогласования Е 
соответственно на коэффициенты усиления С2 и С1. Составляющую $1 сигнала коррекции 
ускорения S7 формируют интегрированием сигнала E и умножением результата на 
коэффициент усиления G3. Сигнал $1 обеспечивает компенсацию постоянной или 
медленно изменяющейся погрешности в сигнале AZ (т.е. погрешности типа "дрейф нуля 
акселерометра"). 

Анализ структурной схемы прототипа показывает, что он обладает следующим 
недостатком. Для того, чтобы на выходе второго интегратора (блока 225) был 
сформирован сигнал барометрической высоты (PZ), необходимо, чтобы на его вход 
поступал сигнал производной по времени от той же барометрической высоты E 7) 


dt 
Однако, как следует из схемы, вместо сигнала ар) на вход второго интегратора 


dt 

nocrynaer сигнал вертикальной скорости VZ, представляющий собой производную по 
времени от геометрической высоты (как результат интегрирования вертикального 
ускорения AZ в первом интеграторе). Таким образом, известный способ дает ошибку, если 
расхождение между геометрической и барометрической высотой становится существенным 
(например, при влиянии метеоусловий). 

Задачей изобретения является повышение точности измерения барометрической 
высоты и вертикальной скорости при полете ЛА в условиях нестандартной атмосферы. 
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Для решения задачи предлагается способ определения барометрической высоты M 
вертикальной скорости, в котором измеряют вертикальное ускорение и с помощью датчика 
статического давления приборную барометрическую высоту, значение вертикальной 
скорости получают путем интегрирования разности вертикального ускорения и поправки 
ускорения, значение барометрической высоты получают путем интегрирования разности 
вертикальной скорости и поправки скорости, разность между полученным таким образом 
значением барометрической высоты и значением приборной барометрической высоты 
используют для вычисления упомянутых поправок ускорения и скорости, отличающийся 
тем, что при определении упомянутой барометрической высоты перед нахождением 
разности вертикальной скорости и поправки скорости вертикальную скорость умножают на 
безразмерный множитель, представляющий собой отношение температуры воздуха по 
стандартной атмосфере к фактической температуре воздуха, измеряемой датчиком 
температуры. 

Изобретение поясняется фигурами. 

На фиг.1 изображена структурная схема алгоритма, лежащего в основе известного 
способа. 

На фиг.2 изображена структурная схема алгоритма, используемого при реализации 
предлагаемого способа. 

На фиг.3 показаны зависимости погрешностей измерения барометрической высоты 
предлагаемым и известным способами для фрагмента полета. 

На фиг.4 показаны зависимости погрешностей измерения вертикальной скорости 
предлагаемым и известным способами для фрагмента полета. 

Для того чтобы устранить указанное несоответствие между барометрической и 
геометрической высотами, сигнал VZ, поступающий на вход второго интегратора, должен 
быть предварительно преобразован в сигнал f діру) . Необходимая функциональная 


dt 
зависимость для преобразования У? B f агро); может быть найдена с использованием 


dt 


следующих исходных соотношений: 
- уравнение статики атмосферы в дифференциальном виде 
ЧЕ Р 


где P, Тн - давление и температура воздуха (по шкале Кельвина), 

К - геометрическая высота, 

R-287,05287 Дж.кг".К`1 - удельная газовая постоянная, 

д - ускорение свободного падения (зависит от широты и высоты ЛА); 

- дифференциальное уравнение, используемое в авиации для формирования 
барометрической высоты полета ЛА в следующем виде 
често @ 
dH, ЕТ 

где Hg - барометрическая высота полета ЛА, 

9с=9,80665 м/с? - стандартное ускорение свободного падения, 

Tcr - температура воздуха по стандартной атмосфере. 

Уравнение (1) следует из формул (1) и (2), приведенных в приложении к ГОСТ 4401-81, 
стр.171, раздел "Уравнения статики атмосферы и состояния идеального газа", уравнение 
(2) для формирования барометрической высоты полета ЛА (Hg) получается из уравнения 
(1) заменой параметра g на gc и параметра Ty на Ter. 

Температуру воздуха по стандартной атмосфере (Тст) определяют в зависимости OT 
барометрической высоты (см. там же стр.174, 175, раздел "Температура, вертикальный 
градиент температуры") следующим образом (до высоты 20000 м): 

T, = 288,15- 0, 0065H,, при - 2000 < H, < 11000 м, (3) 
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Ta = 20665 при 11000 < Н, < 20000 м. (4) 


Из уравнений (1) и (2) путем исключения давления получается следующее уравнение, 
связывающее производные по времени от геометрической и барометрической высот: 


ан т. а dh 
dt — TH g, dt 
В уравнении (5) отношение ад в контексте рассмотренной задачи можно заменить 


Cr 
единицей ввиду того, что его отличие OT единицы является 


малым до высоты 20000 м на всех географических широтах). С 
(— - Ц < 0005 
dr 


учетом сказанного, а также принимая BO внимание TO обстоятельство, что производная по 


времени от геометрической высоты фан у представляет собой вертикальную скорость ЛА 
dt 


(Му B принятых B нашей стране обозначениях, т.к. ось, направленная вертикально вверх, 
обозначается у), уравнение (5) преобразуется к следующему виду: 
ЯН T, 
в дує. 49 
dt qm 


Таким образом, для преобразования Vy B dH, сигнал Му должен бьть умножен на 


dt 
представляющий собой отношение температуры воздуха по 


безразмерный множитель т.. , 


T, 


стандартной атмосфере к температуре воздуха Ha текущей высоте полета, измеряемой 
датчиком температуры. Отсутствие такого множителя приводит к увеличению 
погрешностей в сигналах барометрической высоты и вертикальной скорости при полете в 
условиях нестандартной атмосферы, т.е. в условиях, когда температура воздуха Тн на 
текущей высоте полета будет значительно отличаться от температуры воздуха по 
стандартной атмосфере Тст. 

Предлагаемый способ заключается в следующем (см. фиг.2). 

Измеряют фактическую температуру воздуха Ty, вертикальное ускорение Ау и 
приборную барометрическую высоту нт ‚ причем нх измеряют с помощью датчика 


статического давления (см. приложение к ГОСТ 3295-73, стр.74, формулы 1.1.1.-1.2.2. 
для вычисления геопотенциальной высоты по давлению). Из величины вертикального 
ускорения Ау вычитают (в блоке 1) поправку ускорения $7 и путем интегрирования этой 


= 


разности (в блоке 2) получают значение вертикальной скорости ТЕ которое затем 
Y 


умножают (масштабируют) на множитель т,, (формируемый B блоке З). Из результата 


Ти 


умножения (т.е. из масштабированной вертикальной скорости) вычитают поправку 
вертикальной скорости 55 (в блоке 4) и путем интегрирования разности (в блоке 6) 


получают выходной сигнал (значение) барометрической высоты ја 
Е 


Поправки ускорения 57 и вертикальной скорости 55 находят с использованием 
формируемой в блоке 7 разности высот Е между вычисляемой барометрической высотой 


г и приборной барометрической высотой н” . При этом поправку ускорения 57 находят в 
h 


блоке 8 как сумму умноженной на коэффициент G2 разности высот E и умноженного на 
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коэффициент G3 результата интегрирования B блоке 9 тои же разности высот, в качестве 
поправки вертикальной скорости S5 принимают умноженную на коэффициент G1 ту же 
разность высот Е. 
Температуру воздуха по стандартной атмосфере Тст вычисляют в блоке 5 по формулам 
(3) и (4) с помощью функции [ | ‚ обеспечивающей вычисление Тст. 
b 


Способ может быть реализован как C помощью аналогового, так и цифрового 
устройства. 

Если используется аналоговое устройство, оно должно реализовывать схему, 
изображенную на фиг.2. Как видно из этой схемы, на один из входов блока 3, 
осуществляющего вычисление множителя т,, , подается формируемая блоком 5 функция 


T, 


. Таким образом, хотя до момента определения 
Ћ 


A 


определяемой барометрической высоты В 


A 


значение барометрической BBICOTbI ‚ еще не известно, из-за обратной связи, 


h 
сформированной замкнутой цепочкой блоков 3-4-6-5-3, для определения 


множителя Т,, используется величина определяемой барометрической высоты - 


T, 


Аналогичные обратные связи имеют место при определении поправок скорости и 
ускорения. Поправку скорости находят благодаря обратной связи, образованной цепочкой 
блоков 4-6-7-10-4, а поправку ускорения - обратной связи, образованной цепочкой 
блоков 1-2-3-4-6-7-(11 и 9-12)-8-1. 

Если используется цифровое устройство (компьютер), интегрирование выполняется с 
некоторым временным шагом, начиная с момента, например, взлета летательного 
аппарата. При этом процедуры интегрирования могут быть выполнены на основе 
следующих итерационных формул (эти формулы реализуют схему, изображенную на 


фиг.2), связывающих значения параметров барометрической высоты Е ‚ вертикальной 
Е 


h 


скорости s и необходимой для зтих целей интегральной поправки 
y 


E 


ускорения (обозначаемые далее ^ 


Нв, "Vy, ААУ, 
друга на интервал времени т дискретные моменты 


ААУ ) в равноотстоящие друг от 


^ 


времени киј(_  ,. ): 
НЕ, +, Li Ук+1 А А Урат 


А А Pt А А 
Нь, = Нв, + ски У War + Ay abc B AR yr ij J = ci[Ph, | = н; } (7) 


xi 


2 А A 


У Уха Я М Ук t Рут“ ДАТ Е сар . E на] (8) 


ААу,, = З AY, + c3[ph, . = Hika ў (9) 

где Ріст - прогноз температуры воздуха по стандартной атмосфере, вычисляемый по 
приведенным далее формулам, 

Phg - прогноз барометрической высоты, вычисляемый по приведенной далее формуле, 

t=tk+1-tk - интервал времени. 

Для выполнения цифрового интегрирования по формулам (7), (8) и (9) приближенно 
задают начальные значения параметров на нулевой вычислительной 


Hh "МУ uA AY, 
Hua Mule EEUU 


(при к=1,2,3,...) выполняют с использованием значений этих параметров на предыдущей K- 
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й итерации, измерений температуры наружного воздуха T, , , вертикального ускорения 
ETI 


E 


Ау ки приборной барометрической высоты X а также с использованием прогнозов 


Бун 


температуры воздуха по стандартной атмосфере Pt и барометрической 
"Та 


высоты Ph предварительно вычисленных по следующим формулам: 


Бучі! 


, при 


28815 - 040065 На, - 200 < Hg < 11000м, 


- при : 
а РА ПРИ ocn e «C пала 


я 


Pt. = P Teresa У Av /2-AA VT /2 
в... = HET Е мит + Аут / 2 - YT F . 
Не 


Т.о., при использовании цифрового устройства осуществляется достаточно точная для 
целей определения множителя т,. оценка прогнозного значения барометрической 


Ty 


высоты, что практически равноценно описанной выше обратной связи B аналоговом 
устройстве. 

Сравнение известного и предлагаемого способов в условиях нестандартной атмосферы 
отражено на графиках, приведенных на фиг.3 и 4, следующими параметрами: 

D H1, D Vy1 - погрешности измерения барометрической высоты и вертикальной 
скорости предлагаемым способом, 

D H2, D Vy2 - погрешности измерения барометрической высоты и вертикальной 
скорости способом, реализуемым в прототипе, 

Vy - вертикальная скорость полета ЛА, 

D Тнест - разность температур Тст-Тн, принимаемаа за меру нестандартности 
атмосферы. 

Из фиг.3 и 4 следует, что отсутствие безразмерного множителя Тар В известном 


T 


способе приводит (при относительно небольшой нестандартности атмосферы, 
достигающей примерно 13 градусов) к погрешности измерения вертикальной скорости 74 
м/с, барометрической высоты 712 м (по модулю). При нестандартности атмосферы B 
пределах 20...30 градусов (что достаточно часто имеет место в реальных условиях) 
ожидаемое значение погрешности измерения вертикальной скорости на аналогичных 
режимах полета ЛА достигнет до 10...15 м/с. То есть погрешности определения 
упомянутых параметров уменьшаются (по сравнению с приведенными погрешностями для 
прототипа) до значений 6 M по высоте (параметр D H1 на фиг.3) и 0,9 м/с no 
вертикальной скорости (параметр D Vyl Ha фиг.4) на том же режиме полета ЛА. 

Таким образом, предлагаемый способ позволяет повысить точность определения 
барометрической высоты и вертикальной скорости в условиях полета ЛА в нестандартной 
атмосфере. 


Формула изобретения 
Способ определения барометрической высоты и вертикальной скорости, при котором 

измеряют вертикальное ускорение и с помощью датчика статического давления приборную 
барометрическую высоту, значение вертикальной скорости получают путем интегрирования 
разности вертикального ускорения и поправки ускорения, значение барометрической 
высоты получают путем интегрирования разности вертикальной скорости и поправки 
скорости, разность между полученным таким образом значением барометрической высоты 
и значением приборной барометрической высоты используют для вычисления упомянутых 
поправок ускорения и скорости, отличающийся тем, что при определении упомянутой 
барометрической высоты перед нахождением разности вертикальной скорости и поправки 


скорости вертикальную скорость умножают на безразмерный множитель, представляющий 
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собой отношение температуры воздуха по стандартной атмосфере к фактической 
температуре воздуха, измеряемой датчиком температуры. 
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‚ Му 


D H1,D H2 


44960 44970 44980 44990 45000 45010 45020 45030 45040 
Время, c 
Фиг. 3 
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— D му, м/с 
——D Vy2, міс 


——0 Тнест, Град | 


44960 44970 44980 44990 45000 45010 45020 45030 45040 
Время, с 
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